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ABSTRACT 

We successfully detected a molecular outflow with a scale of 370 - 450 pc in the central region of the 
starburst galaxy NGC 3628 through deep CO(1-0) observations by using the Nobeyama Millimeter 
Array (NMA). The mass of the outflowing molecular gas is ~ 2.8 x 10 7 M Q , and the outflow velocity 
is ~ 90±I0 km s _1 . The expansion timescale of the outflow is 3.3 - 6.8 Myr, and the molecular gas 
mass flow rate is 4.1 - 8.5 M Q yr _1 . It requires mechanical energy of (1.8 - 2.8) x 10 54 erg to create 
this sub-kpc scale molecular outflow. In order to understand the evolution of the molecular outflow, 
we compare the physical properties between the molecular outflow observed from our NMA CO(1-0) 
data and the plasma gas from the soft X-ray emission of the Chandra X-ray Observatory (CXO) 
archival data. We found that the distribution between the molecular outflow and the strong plasma 
outflow seems to be in a similar region. In this region, the ram pressure and the thermal pressure 
of the plasma outflow are 10~( 8 ~ 10 ) dyne cm -2 , and the thermal pressure of molecular outflow is 
10~' 11-13 ) dyne cm" 2 . This implies the molecular outflow is still expanding outward. The molecular 
gas consumption timescale is estimated as 17 - 27 Myr, and the total starburst timescale is 20 - 
34 Myr. The evolutionary parameter is 0.11 - 0.25, suggesting that the starburst activity in NGC 
3628 is still in a young stage. 

Subject headings: ISM: bubbles — ISM: jets and outflows — galaxies: individual (NGC 3628) - 
galaxies: ISM — galaxies: starburst 



1. INTRODUCTION 

Molecular gas consumption, mainly dominated by star 
formation and losing gas through outflow, plays an im- 
portant role on star formation. To understand the molec- 
ular gas consumption in galaxies help us to study galaxies 
evolution. 

Galaxies with violent star formation, especially in 
central nuclear region, are called starburst galaxies. 
They have star formation rates (SFRs) several to hun- 
dreds times higher than normal spiral galaxies. The 
strong star formation activities of starburst galaxies 
produce a large number of young and massive OB 
stars in a very short time. These massive stars fi- 
nally end up their lifetime with supernova explosions. 
Therefore, vast amount of stellar winds from mas- 
sive stars and supernova explosions generate huge en- 
ergy and high pressure to create high velocity galac- 
tic winds, which sweep up the surrounding interstel- 
lar medium (ISM) and create galactic scale outflows 
or superbubbles. These scenarios h ave been described 
by s e veral analytical models ( e.g.. iChevalier fc Cleggl 
[19851 IMcCrav fc Kafatosl [j~987l lYokoo et al.1 119931). nu- 



plasma outflows, such as Ha emission observat ions (e.g., 
lArmus et all 119901: iLehnert fc Heckmanl 1 19961) and soft 
X-ray observations (e.g.. iStrickland et al.ll2004f l. How- 
ever, the observations on molecular outflows or super- 
bubbles are very rare in the past. Lack of this study 
will cause under-estimation on the molecular gas con- 
sumption in galaxies, and have a bias to understand the 
galaxies evolution in our universe. 

The reasons for the rarity of molecular gas detections 
were mostly due to their diffuse/extended nature and 
poor instrumental sensitivities. Recent improvements of 
various instruments provide a better chance to study the 
molecular outflows and superbubbles in galaxies. The 
recent techniques to detect molecular outflows and su- 
perbubbles in galaxies can be classified as three types: 
(1) High- velocity CO wings to detect mole cular out- 
flows, e.g., NGC 3256 (ISakamoto et al.H2006blL M rk 231 
(jFeruglio et al.l [2010h. M 5 1 ([Matsushita et all l2007h . 
NGC 1266 (lAlatalo et all |2011|), and local ULIRGs 
(|Chung et al.ll201lD . (2) P-cvgni profile to detect molec - 
Arp 220 (jSakamoto et alll2009f >. 
MM). 



ular superbub bles, e.g., 
and Mrk 231 (Fis her et al.ll2010D . (3) Direct imaging at 



merical simulations (e.g.. iTomisaka fc Ikeuchil 119881: 
Cooper et al.l I2008D. and observation al results (e.g., 
Heckman et al.l Il990t iCecil et al.l 120021 ). Previous ob- 



servations mainly focused on the warm ionized and hot 



millimeter waveband for CO molecular out flows and su- 
perbubbles, e.g., molecu lar outflow in M82 (jNakai et al.1 
ll987tlWalter et al Il2002h. molecular s uperbubble in M82 
(jWeiss et al.l 119991: iMatsushita et al.l 120001 I2005D. and 
molecular superbubbles in NGC 253 (jSakamoto et al.l 
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2006). The first and second techniques are useful to 
detect outflows and superbubbles from face-on galaxies. 
However, it is difficult to distinguish the molecular gas of 
outflows/superbubbles from that of galactic disk. There- 
fore, we cannot have an accurate measurement on molec- 
ular gas mass losing. To solve this problem, we need to 
observe edge-on galaxies by using the third technique. 
Although this technique is time consuming, we can di- 
rectly observe the morphology of molecular outflows and 
superbubbles, and therefore can have an accurate mea- 
surement on molecular gas mass losing. 

So far, only a few molecular outflows and superbubbles 
in galaxies have been directly imaged by using the third 
technique. More observations to directly image molec- 
ular outflows and superbubbles are necessary in order 
to understand the general properties of these structures 
and their influence on starburst activities. We have been 
conducting detail studies of one molecular outflow and 
two superbubbles t oward the nearb y edge-on starburst 
galaxy NGC 2146 (|Tsai et all [2009ft . Here, we provide 
another sample, NGC 3628. 

Starburst galaxy NGC 3628 is a member of the in- 
teracting Leo Triplet (Arp 317), including NGC 3627, 
NGC 3623, and NGC 3628. NGC 3628 is the north 
source in the Leo Triplet, NGC 3627 is the south one, 
while NGC 3623 is the southwest one. Two signatures 
of tidal interaction in Leo Triplet are observed in several 
wavebands: A plume extending ~ 100 kpc from the east- 
ern edge of NGC 3628 is detected in op tical, far-infrared, 
and HI emission ()Chromev et al.|[T998ft ; a bridge between 
NGC 3628 and NGC 3627 is detected in HI emission 
dCole et all 11998ft . These structures have been modeled 
bv lRotsI (11978ft" and were explained as a tidal interaction 
between NGC 3628 and NGC 3627 about 800 Myr ago. 

NGC 362 8 is a nearby (D = 7 7 Mpc; \" = 37 pc; 
PTulrvl 11988ft . edge-on [i = 87°; iTulTvl 11988ft. a nd IR 
luminous (log Lir/L© = 10.25; ISanders et all 12003ft 
galaxy. The optical image of NG C 3628 looks like a 
boxy shape dChromev et all 119981 see Figure QJ. Ra- 
dio obsexvations~JCondo^ll98^show that NGC 3628 
has a circumnuclear starburst. Besides, a large-scale 
galactic wind has been detected in NGC 3 628. Obser- 
vations with ROSAT (iDahlem et all 11996ft and Chan- 
dra X -ray Observatory (CXO) (|Strickland et all I2001L 
2004) show that NGC 3628 has an asymmetric plasma 
bipolar outflow with ~ 7 - 10 kpc scale in both north 
and south parts; optical Ha observations also show 
that NG C 3628 has a ~ 1 kpc -scale warm plasma gas 
outflow ( Stricklan d et all [2004) . It is natural to ex- 
pect t he existence o f molec ular outflows or superbub- 
bles. Ilrwin fc Sofui (| 1996ft claimed that their CO(l- 
0) Nobeyama Millimeter Array (NMA) observations de- 
tected four expanding molecular superbubbles, which are 
associated with the low velocity gradient ridge outside of 
the nuclear disk. However, no molecular outflow has been 
detected. Our new CO(1-0) NMA observation detect a 
sub-kpc scale molecular outflow for the first time. We 
will discuss its properties, its impact on NGC 3628, and 
its evolution by comparing with the CXO archival data. 



2. OBSERVATIONS AND ARCHIVE DATA 
2.1. NMA CO (1-0) Observations 



Figure 1. DSS Optical Image of NGC 3628. The box size is 
2' X 2', the same size of NMA CO(1-0) moment maps showing in 
Figure [3] 

We have used the NMA to observe the CO(1-0) emis- 
sion of NGC 3628 in the central ~ V region. The obser- 
vations were made during 2000 December to 2002 Febru- 
ary with 3 different configurations of six 10-meter anten- 
nas. The total on-source time was ~ 40 hours. The 
phase tracking center was a(J2000) = ll h 20 m 17:011 
and (5(J2Q OO) = 13°35'20'/06 4. We used tunerless SIS 
receivers (jSunada et all 11994ft . and observed the CO(l- 
0) line in the upper side band. Double side band sys- 
tem temperature was about 200 - 300 K for most of 
the observations. The backend used was the XF-type 
spectro-correlator Ul tra Wide Band Correlator (UWBC; 
lOkumura et alll2000ft . configured to have 512 MHz band- 
width with 256 channels (i.e., a 2 MHz channel width or 
5.2 km s _1 velocity resolution). We observed 3C 273 and 
3C 279 as bandpass calibrators. The phase calibrator was 
1055+018, and amplitude calibrators were 1055+018, 
3C 84, and 3C 345. The flux scale of 1055+018 in 2000 - 
2001 was determined by comparison with Mars. That in 
2001 - 2002 was determined by comparison with 3C 84 
and 3C 345, which were determined by comparison with 
Mars, Uranus, and Neptune. The uncertainties in the 
absolute flux scales are estimated as 5 - 6% in 2000 - 
2001 and 10 - 12% in 2001 - 2002. 

The data were calibrated by using the NRO soft - 
ware package "UVPROC II" (|Tsutsumi et all [l997ft . 
and were CLEANed by using standard procedures im- 
plemented in the NRAO Astronomical Image Process- 
ing System (AIPS). The maps were made with nat- 
ural weighting, with a final synthesized beam size of 
3'.'01 x 2'.'36 (147 pc x 115 pc), which has a higher reso- 
lution than the past C O (1-0) observations ( 3"93 x 2"79 
or 128 pc x 123 pc; Ilrwin fc Sofui ll~996ft . The po- 
sition angle of the beam is 143°. The noise level 
is 9.86 mJy beam -1 with a velocity resolution of 
5.2 km s" 1 , b etter than the previous observation of 
Ilrwin fc Sofud (|1996ft . 65 mJy beam -1 with a velocity 
resolution of 13 km s" 1 . 

2.2. CXO Archive Data 
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We obtained X-ray data from the CXO archive, which 
were originally observed by Strickland et al. in 2000 De- 
cember 2 with 60 ks total exposure time. It was ob- 
served with the Advanced CCD Imaging Spectrometer 
(ACIS). The nucleus of NGC 3628 is placed in the ACIS- 
S3 chip, which is a back illuminated CCD on the spectro- 
scopic array. The back illuminated CCD is more sensitive 
to low-ener gy X-ray photons than the front-illumi nated 
ACIS chips tiStrickland et al.ll200l llnui et al1l2005l) . and 
therefore is powerful to observe diffuse X-ray emission 
in low energies, such as the X-ray superwinds. The dif- 
fuse gas images of NG C 3628 have been published in 
IStrickland et al.l ((2004), but the detail properties have 
not yet been discussed. In order to understand the prop- 
erties of the diffuse gas, such as density, mass, energy, 
and pressure, we reprocessed the data, extracted the 
spectra, and performed model fittings. 

The data were reprocessed by using the Chandra Inter- 
active Analysis of Observations software package (CIAO) 
version 4.2, and the Chandra Calibration Database 
(CALDB) version 4.2.2, released on 19 April 2010. 
We used the High Energy Analysis software (HEAsoft) 
version 6.5, and the X-Ray Spectral Fitting Package 
(XSPEC) version 12.4 for further analysis. The back- 
ground is chosen from a source-free circular region with 
a radius of l'.3 centered at the position a(J2000) = 
ll h 20 m 04 i !86 and <5(J2000) = 13°30'50"49, ~ 5' away 
from the galactic center. The background level was 0.04 
counts s arcmin -2 . To avoid the bias of choosing back- 
ground, we double checked with blank-sky subtraction. 
We used CIAO script DEFLARE and chose sigma as the 
flare-cleaning method. After clipping data with count 
rates less than 3er, the exposure time with good time 
intervals turned to be 56,319 s. 

We only consider the data with energy between 0.3 keV 
and 7.0 keV. This is because the data with energy above 
7.0 keV is dominated by X-ray background, and the data 
quality below 0.3 keV has large uncertainty due to poor 
calibration. We grouped the energy bins of the spec- 
tra from the selected region so that each bin would con- 
tain at least 20 counts to allow \ 2 statistics for spec- 
tral fitting. In order to analyze the diffuse and extended 
X-ray emission, we used two CIAO threads, WAVDE- 
TECT and CELLDETECT, to detect point sources and 
removed them from the X-ray data. After removing de- 
tected point sources from the X-ray data, we used the 
CIAO thread DMFILTH to interpolate the emission at 
point source regions from their surrounding backgrounds. 

3. RESULTS 
3.1. NMA CO(1-0) data 
3.1.1. Distribution 

We summed every 4 channels (= 20.8 km s -1 ) from the 
original channel maps to reduce the rms noise level. The 
CO(1-0) channel maps are shown in Figure[5J The noise 
level is 4.93 mJy beam -1 (68.8 mK). Comparing to the 
noise level of llrwin fc Sofuel (fl99ll . 65 mJy beam with 
a velocity resolution of 13 km s , which corresponds to 
52 mJy beam -1 with a velocity resolut ion of 20.8 km s -1 , 
the emission thatE rwm fc Sof uc (1996!) detected is - 10ct 
and the second contour in our Figure O This suggests 
th at we could detect 10 times weaker emission than that 
of|l rwin fe Sofuel (|1996l) . The emission in each channel 



shows that most of the gas is distributed along the galac- 
tic disk, and the structure in either one or both edges of 
the galactic disk is fragmented. Besides, some extended 
diffuse features appear in the north of the galactic center 
in several channels. 

Figure [3^ is our NMA CO(l-O) integrated intensity 
map (moment zero map). Most of the CO emission is 
distributed along the galactic disk with a position an- 
gle of 104°. The distribution of molecular gas in the 
galactic disk is not symmetric with respect to the major 
axis. The amount of the molecular gas above the major 
axis is more than that below the major axis. Besides, 
both edges of the molecular disk are distorted. The east- 
ern edge of the disk shows a boxy-shape structure, and 
the western edge of the disk is fragmented. Both edges 
looks similar to the edges of the optical image (Figure [lj 
although the scale is different. An extended emission lo- 
cated from the center with an offset of (AR.A., ADec.) 
~ (0", +15") labeled as OF is the real emission because 
they also appear in the same position by mapping with 
uniform weighting. However, those centers with position 
offsets of - (0", +25") and - (+7", -20") are considered 
as sidelobes because the positions of these features are ei- 
ther disappeared or shifted along with the beam patterns 
by mapping with uniform weighting or by choosing dif- 
ferent UV-ranges. The existence of the strong sidelobes 
is due to the UV-coverage of a low-declination source, 
which generates strong sidelobes along the north-south 
direction. 

3.1.2. Velocity Features 

Figure [5Jd is the intensity weighted mean velocity map 
(moment one map). From this map, we can see that 
there is a rigid-body rotation in the galactic center and 
a flat rotation in larger radii in the eastern part of the 
disk. The velocity in the northwestern part of the disk 
seems to be a flat rotation, but the southwestern part 
of the disk does not show a clear flat rotation feature. 
Besides, the velocity in the western part of the disk is 
twisted and entangled. In order to further discuss the 
detail of these features, we first rotate the moment zero 
map until the major axis of the galactic disk is the same 
with the horizontal direction, namely rotating clockwise 
for 14°, as shown in Figure HJ and then we check the 
position-velocity (hereafter p — v) diagrams parallel or 
perpendicular to the major axis. 

Figure [5] shows the major-axis p — v diagram by averag- 
ing the p — v diagrams within 2" (between +1" and —1" 
offsets) from the major-axis. From Figure [SJ we found 
that the galactic disk consists of several structures: The 
curve between the position +60" and —60" indicates a 
rotation curve along the major axis of the galactic disk, 
namely, a rigid-body rotation between the positions +10" 
and —10", with a flat rotation beyond these positions; 
this suggests the existence of a large-scale molecular gas 
disk. The solid line between the positions +4" and —4" 
indicates a rigid-body rotation, suggesting the existence 
of an inner molecular disk. The ellipse between the posi- 
tions +40" and —40" indicates the velocity feature more 
or less symmetric. Th is feature is similar to the p — v di- 
agram of our Galaxy (|Englmaier fc Gerhard! 1 19991). sug- 
gesting the existence of a molecula r bar ( Binnev et al.l 
Il99lt iGarcia-Burillo fc Guelin|[l995] ) in this galaxy. The 
solid line between the position +40" and —40" will be 
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Figure 2. NMA CO(l-O) channel maps of NGC 3628. The red cross is the phase tracking center of our observation, «(J2000) = 
ll h 20 m 17?011 and <5(J2000) = 13°35'2CK'064. The synthesize beam size is 3'.'01 X 2'.'36 (112.4 pc X 88.1 pc), which is shown in the lower 
left corner of the first panel. The contour levels are 4, 10, 20, 40, 60, 100, and 130<r, where la is 4.93 mjy beam -1 . The LSR velocity, in 
units of km s — 1 , is showing at the top left corner of each channel map. 
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Figure 3. The NMA CO(l-O) integrated intensity (moment zero) and intensity weighted mean velocity field (moment one) maps of NGC 
3628. The cross is the phase tracking center and the synthesized beam is shown in the bottom-left corner of each map. The central position 
and the synthesized beam size are the same as in Figure [2] The red dashed line indicates the location of the galactic disk, and the red 
curve indicates the molecular outflow OF. (a) Moment map. The contour levels are 1, 3, 5, 10, 20, 50, 100, 150, and 180<r, where lcr is 
648.3 mjy beam -1 km s' 1 (= 9.05 K km s" 1 ). (b) Moment 1 map. The contour levels are from 640, 660, 680, . . ., and 1,040 km s , 
increasing with 20 km s —1 . 
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Figure 4. NMA CO(l-O) moment zero map after clockwise ro- 
tation of 14°. The contour levels are the same as in figure [3^. The 
two horizontal dashed lines, offset by +12" and —12" from the 
major axis, specify the range of major-axis p — v diagrams shown 
in Figure|S] The two vertical dashed lines, offset by +1" and —11" 
from the minor axis specify the range of OF. Its minor-axis p — v 
diagram is shown in Figure [7] 



to ±40", so that this structure may be a rotating disk 
which is located at a large radius. Here we call it an 
outer molecular disk. 

In order to see the velocity feature of the extended 
structure OF in the north of the galactic disk, we av- 
eraged the p — v diagrams along the minor axis in the 
range between the position offsets +1" and —II" on the 
major axis (see Figure [4] for this range), and the aver- 
aged p — v diagram is shown in Figure [7] The strong 
emission is distributed along the position offset of ~ 0", 
indicating the location of the galactic disk. Northern 
part of the galactic disk (positive side of the position off- 
set) obviously shows an extended diffuse emission, which 
corresponds to the extended structure OF in Figure [3^ 
and Figure |4j The velocity range of this diffuse emission 
is between ~ —110 — h-70 km s _1 . On the other hand, 
southern part of the disk (negative side of the position 
offset) shows no extended emission. 
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Figure 5. The p — v diagram along the major axis by averag- 
ing the data between the positions —1" and +1" from the major 
axis. The zero velocity corresponds to the systemic velocity of 
840 km s" 1 . The contour levels are 2, 3, 5, 10, 15, 20, and 50cr, 
where la is 9.09 mjy beam -1 . The solid curve between +60" and 
—60" indicates the rotation curve of the galactic disk. The solid 
line between two dashed vertical lines at +4" and —4" indicates 
the steep velocity gradient of the inner molecular disk. The solid 
line between two dashed vertical lines at +40" and —40" indicates 
the low velocity gradient of the outer molecular disk. The solid 
ellipse indicates the velocity feature of the molecular gas bar. 

discussed in the following paragraph. 

Figure [B] shows averaged p — v diagrams parallel to the 
major axis with various offsets. These p — v diagrams 
are obtained by averaging every 2" in the range between 
the position offsets +12" and —12" from the major axis 
(see Figured] for this range). The red solid curve, lines, 
and ellipse are the same as those in Figure [Sj We found 
that the structure of the molecular disk is not symmetric 
in the north-south direction. The inner disk can only be 
clearly seen in the position at the range between the offset 
+1" and —1". The velocity feature of the bar structure 
can be seen in the position at the range between the offset 
+7" and —5". In addition to these velocity features, we 
also found a rigid-rotation feature in the position at the 
range between the offsets —7" and —9". The emission is 
weak while the slope is shallow, and the features extend 



3.1.3. Mass 

The total molecular gas mass can be estimated from 



the CO flux. The total CO flux, S, 



> J CO(1-0)) 



measured 



from our moment zero intensity map (Figure [3Ji) is ~ 
8.01 x 10 3 Jy km s" 1 (1.12 x 10 5 K km s" 1 pc 2 from the 
central w 80" x 20" of the galaxy) The total flux is four 
times larger than that measured bv llrwin" fc Sofud (fl996h 
(~ 2.05 x 10 3 Jy km s _1 ) from the same region of the 
galaxy. This is because our high sensitivi ty NMA obser- 
vation provides a better sensitivity than (irwin fc Sofud 
( 1996) (as mentioned in Sec. 13.1.11) so that we can detect 
much weaker emission. The total molecular gas mass de- 
rived from our CO data, in units of M , can be calculated 
as 



M H2 =1.2 x 10 4 x D 2 

x Sco(i-o) x 3 _ 0xl()20 , 
M- B8 = 1.36x M H3 , 



(1) 

(2) 

where D is the distance in units of Mpc, Xqo is the CO- 
to-H2 conversion factor in units of cm -2 K (km s -1 ) -1 , 
and the factor 1.36 is a coefficient convert from molecular 
hydrogen gas mass to total gas mass, including helium 
([Sakamoto et al.lfl995[ ). Here, we used the value of Xco 
= 1.4 x 10 20 cm" 2 K km s" 1 ([Matsushita et al.ll2000h . 
assuming the condition of the molecular gas is similar to 



M„ 



1.8 x 



that of M82. Thus the total gas mass, 
10 9 M Q . 

The dynamical mass, Mdyn within the radius of r, can 
be derived from the rotational curve. Figure[5]shows that 
the maximum rotational velocity, v rot ~ 220 km s _1 , 
appears at the location of radius ~ 40" (1.5 kpc). We 
assume a circular rotation for the galactic disk, v 2 ot — 
GMdy n /r. Therefore, the dynamical mass within 1.5 kpc 
radius can be estimated as ~ 2 x 10 10 M Q . This indicates 
that the gas mass to dynamical mass ratio of 10%, which 
is the typical value for the central regions of galaxies 
(c.f., from 4% for Sa g alaxies to 25% for Scd galaxies; 
lYoung fc ScoviililT99lK . 



3.2. CXO X-ray Archive Data 
3.2.1. Image of the Diffuse Emission 
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Figure 6. Averaged p — v diagrams parallel to the major axis in the range between the position offset of —12" and +12", which are 
marked as two dashed horizontal lines in figure [4] Each p — v diagram is averaged every 2" with an interval of 2" . The number in each 
p — v diagram indicates the offset from major-axis in unit of arcsec. The solid lines, curves, and ellipses are the same as those shown in 
Figure[5] The zero velocity corresponds to the systemic velocity of 840 km s — . The contour levels are 3, 5, 7, 10, 20, 30, 40, and 50cr, 
where lcr is 9.09 mjy beam -1 . 
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Figure 7. The p — v diagram along the minor axis averaged over 
the range between the position offsets +1" and —11" from the mi- 
nor axis, which is between two dotted vertical lines in figure[4] The 
zero velocity corresponds to the systemic velocity of 840 km s -1 . 
The contour levels are 2, 4, 10, 20, 30, 50, and 70cr, where la is 
5.23 mjy beam -1 . 

In order to compare the CXO X-ray data with the 
NMA CO(1-0) data in different scale, we made two types 
of smoothed images. 

To image the large scale diffuse structure, we used the 
CIAO tool ACONVOLVE to convolve the point-source- 
removed X-ray image with a Gaussian that is smoothed 
by 25" in 2 energy bands, soft (0.3 - 2.0 keV) and hard 
(2.0 - 7.0 keV) bands. We also tried to smooth with 
20", but the diffuse emission cannot be displayed clearly. 
The smoothed soft and hard X-ray images are shown in 
Figure |5k. and Figure [Sb, respectively. These two images 
are ve ry similar to those published by [Stricklan d et all 
(|2001h . suggesting that our imaging is consistent with 
theirs. In soft-band X-ray emission (Figure Ek), the 
strongest intensity is very close to the galactic center, 
but shifted a little bit toward the north. An asymmet- 
ric plasma outflow located above and below the galactic 
center is clearly seen. The northern plasma outflow has a 
stronger intensity and wider distribution than the south- 
ern outflow, while the southern outflow extends longer 
distance than the northern one. The diffuse emission 
spreads over ~ 7 - 10 kpc scale in both north and south 
of the galactic disk. In hard-band X-ray emission (in 
Figure [SJj) , the strongest intensity is at the galactic cen- 
ter. The emission is roughly concentrated in the galactic 
center, and we do not see clear evidence of the outflow 
feature. The overlaid images between the CXO X-ray 
data and the NMA CO (1-0) data are shown in Figure^ 
and Figure [5b. 

To image the small scale diffuse structure, we used 
AIPS task CONVL to convolve the X-ray data to the 
same beam size as (3"01 x 2"36) as the NMA CO(1-0) 
data in 2 energy bands. The high resolution soft and hard 
X-ray images are shown in Figure [SJ; and Figure IHJi, re- 
spectively. In the higher resolution soft-band X-ray emis- 
sion map (Figure®:), the distribution of the plasma out- 
flow near the galactic center is roughly located in similar 



region as the extended molecular gas feature OF, while 
the peak intensity of the X-ray emission is located a little 
bit toward the north. In the higher resolution hard-band 
X-ray emission map (Figure®!), the distribution is con- 
centrated in the galactic center. 

3.2.2. Spectral Analysis of the Diffuse Emission 

In order to compare the physical properties of plasma 
gas to that of molecular gas, we used XSPEC to fit the 
spectra of the X-ray data to obtain the temperature kT, 
and the normalization term Norm, which can be used 
to estimate mean electron density, plasma mass, ther- 
mal pressure, and thermal energy (the details will be 
described in Sec. 13.2.3"]) . 

For this spectral analysis, wc first define spectral fitting 
regions. FigurefTUk and FigurellOb sh ow the VLA FIRST 
HI archive image (jBecker et alJll995l ) overlaid on the soft 
and hard X-ray image, respectively. The HI line emission 
is distributed along the galactic disk. After comparing 
the distribution of X-ray emission, HI emission, and the 
extended molecular gas OF, we specified/defined two fit- 
ting regions, CORE and HALO. For simplicity, we de- 
fined the region CORE as a circular area with a radius 
of ~ 0.'3 (0.67 kpc), which includes the central strong 
X-ray emission in the HI disk and the feature OF (Fig- 
ure^, ®i, and HUk) . Meanwhile, we defined the region 
HALO as a circular area with a radius of ~ 3.'0 (6.72 kpc), 
which includes all diffuse X-ray emission but excludes the 
region CORE. 

The temperature kT, and the normalization term 
Norm, depend on the absorption column density and 
the metal abundance of the galaxy. The absorption col- 
umn density of atomic hydrogen, Ah, can be fitted with 
the photoelectric absorption model, WABS. In CORE, 
the absorption is affected by both the host galaxy and 
the Milky Way. Since the absorption in the host galaxy 
is unknown, we left Ah as a free parameter. In HALO, 
X-ray emission has no counterpart in the HI disk, thus 
we only consider the absorption from the Milky Way. 
We therefore fixed Ah as Milky Way absorption col- 
umn density, 0.02 x 10 22 cm" 2 (jKalberla et al.1 [20051 
iDickev fe Lockmanl[l990l ). 

For the metal abundance, we identify O (~ 0.5 - 

0. 7 keV), Si (~ 1.8 - 1.9 keV), and Fe (~ 0.7 - 1.0 keV) 
in the HALO spectrum (Figure [TTV), and Mg (~ 1.3 - 
1.4 keV), Si (~ 1.8 - 1.9 keV), and Fe (~ 0.7 - 1.0 keV) 
in the CORE spectrum (Figure [Tib) . The existence of 
these lines suggests that the emission originat es from 
optically-thin thermal plasma ([Inui "et~aLl [20051) . Thus 
we can choose a thermal-plasma model from XSPEC, 

1. e., VMEKAL or VAPEC. 

Since the data quality of the HALO spectrum is poor 
at energy lower than 0.4 keV and higher than 1.3 keV, 
we intended to exclude these data. In addition to con- 
sidering the metal abundance of Si at ~ 1.8 - 1.9 keV, 
we include data at higher energy side till 3.0 keV. 
Therefore, the soft-band energy range is 0.4 - 3.0 keV. 
We begin to fit the spectrum with an absorbed single- 
temperature thermal-plasma model. However, the spec- 
trum cannot be fitted well neither by WABS*VMEKAL 
nor by WABS*VAPEC model. Thus we tried to fit 
the spectrum via an absorbed two-temperature thermal- 
plasma model, i.e., WABS* (VMEKAL + VMEKAL) or 
WABS* (VAPEC + VAPEC). We also fix the metal abun- 
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where EI is the emission integral, 

EI = I n e nnfdV, 



Figure 8. Smoothed (by 25") CXO X-ray image after removing point sources. The cross is the phase tracking center of the NMA CO(1-0) 
data at a(J2000) = ll h 20 m 04!863 and 5(J2000) = 13°30'50'.'49. (a) Soft band (0.3 - 2.0 keV). The contour levels are 2, 3, 5, 10, and 20 X 
2.6757 counts s _1 . (b) Hard band (2.0 - 7.0 keV). The contour levels are 2, 3, 5, and 10 X 2.5802 counts s" 1 . 

dances of the higher temperature component to those of 
the lower temperature component. The fitting results 
show that the model WABS* (VMEKAL + VMEKAL) 
fits better than the model WABS*(VAPEC + VAPEC). 
The fitting parameters in the model WABS* (VMEKAL 
+ VMEKAL) for both HALO and CORE regions are 
listed in Table □ 

In HALO, the metal abundances fitted by the model 
are ~ 0.5 Z Q for O, - 2.8 Z Q for Si, and - 0.4 Z Q 
for Fe, where Zq is the solar abundance. Two tempera- 
tures are ~ 0.21 keV and 0.60 keV. In CORE, the metal 
abundances fitted by the model are ~ 5.0 Zq for Mg, 
~ 13.5 Zq for Si, and ~ 0.8 Zq for Fe. Two tempera- 
tures are ~ 0.47 keV and 7.3 keV. We notice that both 
the metal abundance of Mg and Si are extremely high. 
These values are likely not physical. Besides, their large 
uncertainties (1.5 - 25.9 Zq for Mg and 5.3 - 50.8 Zq 
for Si) suggest that the fitted values are not accurate 
enough. One possible reason is that the CCD material 
contaminates the detection of Si. Another possible rea- 
son is that the photon counts, ~ 1000 counts, are too low 
to get a reasonable fitting results for Mg and Si. We fix 
the metal abundances of Mg and Si as unity and re-fit the 
spectrum again. This does not change kT and Norm too 
much. Therefore, we could ignored the influence from 
high abundances of Mg and Si. 

3.2.3. The Properties of Plasma Gas 

From kT and Norm, we could derive mean elec- 
tron density n c , plasma mass M p i m , thermal pressure 
-P p im,thm, and thermal energy -E p im,thm- Tne definition 
of Norm is Q 



(4) 



Da is angular diameter distance, z is redshift, V is the 
volume of plasma, Un is the hydrogen density, and / 
is the volume filling factor. For nearby galaxies, Da is 
simply the distance D. Assuming that the electron and 
hydrogen have the same density, i.e., EI = J n^fdV, we 
can derive the following plasma properties. 
The mean electron number density of plasma is, 

The mass of plasma is, 

M p i m = p p i m V = m p n c V, 

where p p i m is the mass density of plasma, and m p 
proton mass. The thermal pressure of plasma is, 

Pplm 



plm,thm " 



kT = 2n r kT, 



fj,m p 



(5) 

(6) 
is the 

(7) 



where \x is the mean molecular weight. For a fully ionized 
hydrogen, there are two particles for every proton and 
\i = 1/2. The thermal energy of plasma is, 



-E'plm.thm — - fplm,thm^ — 3n c VkT 



(8) 



Norm 



10 



47r(£> A (l + z)f 



EI, 



(3) 



1 T he definition is from XSPEC websitej 

http: / /heasarc. gsfc.nasa.gov/docs/xanadu/xspec/manual/ 



(jCarroll fc OstlidflflQfih . 

We applied z — 0.0028 base on the systemic velocity 
measured from our NMA CO(l-O) data. The derived 
plasma parameters from the two-temperature model in 
HALO and CORE regions are shown in Table [2j 

The electron densities of the lower and the higher tem- 
perature component in CORE are 9.5/ -1 / 2 x 10~ 3 cm~ 3 
and 2.0/ -1 / 2 x 10 -2 cm -3 , respectively, about one order 
of magnitude higher than those in HALO. The plasma 
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Figure 9. The spatial correlation between the NMA CO(1-0) data (contours) and the CXO X-ray data (color scale). The contour levels 
are the same as those in Figure[3] The cross is the phase tracking center of the NMA CO(1-0) data. The large and small red solid circles 
are the regions HALO and CORE, respectively. The red dashed curve is the OF. (a) CO(1-0) contours overlaid on the soft band (0.3 — 
2.0 keV) X-ray image, (b) CO(1-0) contours overlaid on the hard band (2.0 - 7.0 keV) X-ray image, (c) Zoom-in image of the central 
2' X 2' region of (a), (d) Zoom-in image of the central 2' X 2' region of (b). 



masses of the lower and the higher temperature com- 
ponent in CORE are 3.0/ 1 / 2 x 10 5 M and 6.4/ 1 / 2 x 
10 5 Mq, respectively, about one to two orders of magni- 
tude smaller than those in HALO. The plasma thermal 
pressures of the lower and the higher temperature com- 
ponent in CORE are 1.4/" 1 / 2 x 10~ n dyne cm and 
4.8/ -1 / 2 x 10~ 10 dyne cm~ 2 , respectively, about one 
to two orders of magnitude higher than those in HALO. 
The plasma thermal energies of the lower and the higher 
temperature component in CORE are 0.8/ 1 / 2 x 10 54 erg 
and 2.7/ 1 / 2 x 10 55 erg, respectively, about one order of 
magnitude smaller than those in HALO. 

3.3. The Extended Structure above Galactic Disk 



Our NMA CO(l-O) observations detected an extended 
structure OF that appears only in the north of the galac- 
tic disk for the first time (Figure [3^,) . After comparing 
the morphology between the NMA CO (1-0) data and the 
CXO X-ray data (Figure [5k, Figured;, and Figured!), 
we found that the distribution of OF is roughly match- 
ing the northern plasma outflow near the galactic cen- 
ter (CORE), and seems to be the ejection point of the 
northern large-scale plasma outflow. The intensity peak 
of X-ray emission is a little bit toward the north of OF, 
which is similar to NGC 2146 (|Tsai et al.|[2009fl . Com- 
bining the kinematic feature of OF in the p — v diagram 
(Figure [7]) , we conclude that OF is a molecular outflow 
and is part of the large-scale HALO outflow. 
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Figure 10. Smoothed (by 25") CXO X-ray image (color) overlaid on the Hi VLA FIRST data (contour). The cross is the phase tracking 
center of the NMA CO(1-0) data, (a) Hi contours overlaid on the soft-band (0.3 - 2.0 keV) X-ray image, (b) Hi contours overlaid on the 
hard-band (2.0 - 7.0 keV) X-ray image. The contour levels are 2, 3, 5, 7, 10, 20, 50, 100, 200, and 500cr, where lcr is 0.15 mjy beam" 1 . 
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Figure 11. The CXO X-ray data with the fitted XSPEC model. The top plots in each figure are the data (crosses) with the fitted model 
(dot for individual model and histogram for combined model). The bottom plots in each figure are x 2 - (a) Spectrum extracted from the 
HALO region with the fitting model of WABS*(VMEKAL + VMEKAL). (b) Spectrum extracted from the CORE region with the fitting 
model of WABS*(VMEKAL + VMEKAL). The fitting parameters are shown in TablcQ] 



The expansion velocity and the size of OF can be es- 
timated from the averaged minor-axis p — v diagram. In 
Figure[71 the velocity range of OF is about ~ 180 km s _1 
(Sec. 13.1.21) . Since we do not see clear velocity gradient 
in a certain direction, the simplest model for OF is an ex- 
panding outflow that explodes isotropically at one side of 
the disk (see the dashed curve in Figure [7]). The expan- 
sion velocity and radius, which measured from the p — v 
diagram, are i> eX p ~ 90± 10 km s _1 and R a f ~ 10" — 12" 
~ 370 - 450 pc, respectively. 

The molecular mass of OF can be estimated from CO 
flux. The total flux, summed from each channel map, is 
82.0 Jy km s" 1 (= 1.14 x 10 3 K km s" 1 ). Using Eq. Q] 
and Eq. [2J assuming the same conversion factor, we can 
derive the gas mass of OF, M oi , as 2.8 x 10 7 M Q . The 
mass is only a few percent of the total molecular gas 
detected in our data. 



The mechanical energy is -E f,mcch — 2 "" 1! ''<>:p 



±M of i4 p = 

(1.8 - 2.8) x 10 54 erg. After considering the mechan- 
ical efficiency 7, ~ 10 - 20 %, the energy transm.it- 
ting from supernova explos i ons into the su r round i ng ISM 
(iMcCrav k. Kafatosl Il987t iWeaver et all 119771: iLarsonl 
11974J). The energy corresponds to 9,000 - 28,000 super- 
nova explosions. Assuming the expansion velocity keeps 
constant, the expansion timescale is t exp = R i/v CX p — 
3.3 - 6.8 Myr. The average molecular gas mass flow 
rate can be calculated after taking account the expan- 
sion timescale, M f = M f/t e x P ~ 4.1 - 8.5 M Q yr _1 . 

4. DISCUSSION 

4.1. Bar, Disk, or Superbubbles? 

Ilrwin fc Sofud (|1996l ) claimed that they detected four 
molecular superbubbles mostly associated with the low 
velocity gradient ridge outside of the nuclear disk. The 
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locations and the expansion velocities of these four 
molecul ar superbubbles are shown in Figure 9 and Ta- 
ble 2 ofllrwin fc Sofud (I1996D . Two superbubbles, A and 
C in llrwin fc SofuT (fl996) , are located at the same po- 
sition and have the same central velocity but different 
expa nsion velocities . Ano ther two superbubbles, B and 
D in|l run fc Sofud (fl99ll . are located at different posi- 
tions and have different central velocities. 

In order to confirm these four structures are real super- 
bubbles from our data, we used two types of images, the 
p — v diagrams (Figure I12p and the intensity maps (Fig- 
ure [T3j) , to examine these structures. We parametrized 
these structures with four parameters, (An, AS, v c , 
v Bxp ), where Aa and AS specify the location, v c is the 
central velocity of the superbubble, and u CX p is the ex- 
pansion velocity of the superbubble. The big blue dot in 
each p — v diagram of Figure [12] is to mark the position 
and the central velocity of each superbubble. The cross 
in each integrated intensity map of Figure II 31 which is 
integrated along superbubble's expansion velocity, is to 
mark the central location of the superbubble. 

Superbubble A has a coordinate of (—13", 0", 
55 km s -1 , 20 km s' 1 ). In the p — v diagram of Fig- 
ure 112b . we do not find any expanding structure at 
position A. In addition, in the intensity map of Fig- 
ure [13k . which is integrated along the velocity range of 
35 - 75 km s _1 , we do not see any extended/distorted 
structure at position A. Since it is difficult to match su- 
perbubble A at position A both in the p — v diagram 
and the intensity map, it is difficult to conclude struc- 
ture A as a superbubble. Meanwhile, the position A in 
Figure [T2k is close to the red ellipse, which represents 
the p — v diagram of a bar structure. This suggests that 
structure A is more likely part of the bar rather than a 
superbubble. 

Superbubble B has a coordinate of (—33", — 1", 
-62 km s _1 , 33 km s' 1 ). Its p — v diagram is shown 
in Figure [12b . and its intensity map, which is integrated 
along the velocity range of 29 - 95 km s _1 , is shown in 
Figure [T3"b . Unfortunately, we do not find any expanding 
structure at position B in the p — v diagram fFigurc [T"2"b ). 
neither find any shell structure at the position B in the 
intensity map (Figure IT3"b ) . Thus it is also difficult to 
identify structure B as a superbubble. Besides, since po- 
sition B in the p — v diagram is located inside the red 
ellipse, it is also possible to consider structure B as part 
of the bar. On the other hand, the emission distribution 
at position B in the p — v diagram looks different from 
that at its symmetric position, which has a coordinate 
of (33", -1", 62 km s -1 , 33 km s" 1 ). Namely, there is 
almost no emission at position B, but there is emission 
detected at its symmetric position. An explanation for 
this asymmetric is that the molecular gas at the posi- 
tion B has been blown up. This indicates a possibility of 
the existence of a superbubble at position B in the past. 
However, we can not distinguish which scenario is true. 
Therefore, it is also difficult to conclude that structure 
B is a superbubble. 

Superbubble C has a coordinate of (—13", 0", 
55 km s , 150 km s _1 ), which has very similar coor- 
dinate with that of superbubble A. The only difference 
is their expansion velocities. Thus we can examine su- 
perbubble C by using the same p — v diagram with that 



of superbubble A (Figure [T^k). Again, we cannot find 
any emission represents superbubble C. Thus it is still 
difficult to identify structure C as a superbubble. 
Superbubble D has a coordinate of (—20", —5", 

- 126 km s' 1 , 85 km s' 1 ). We find emission detected at 
position D in the p—v diagram of Figure fTSb . In addition, 
in the intensity map of Figure 113b . which is integrated 
the velocity range between —211 and —41 km s _1 , we 
find an extended structure at the position D. We further 
check the p — v diagram (Figure I14[) averaged between 
position offsets —16" and —24" along the minor axis of 
Irwin & Sofue's major axis (i.e., clockwise rotation of 
17°) in order to compare the two observations. We de- 
pict their superbubble D in Figure [TJ] by a thin cross 
and a thin dashed ellipse. Interestingly, our deep ob- 
servation detected more diffuse emission. The enclosed 
loop marked as a thick dashed ellipse, which represents 
the structure of a superbubble, is smaller than that of 
Irwin & Sofue's thin dashed loop. We identify the thick 
loop as the superbubble D' with a coordinate of (—20", 
-4", -138 km s" 1 , 50 km s' 1 ). Since a shift of 1" from 
superbubble D is smaller than the angular resolution of 
~ 3", we conclude that the superbubble D and D' are 
identical, but the superbubble D' has a more accurate 
coordinate than the superbubble D. The radius of the 
superbubble D' can be estimated from Figure [T2k and 
Figure [H ~ 5" x 5" (~ 186 pc x 186 pc). 

In conclusion, our results only confirm the existence of 
the superbubble D, superbubble A and C are likely to be 
parts of the bar structure, while superbubble B could be 
either part of the bar structure, or a past superbubble 
whose molecular gas has already been blown up. The 
mis-interpretation is possibly due to the low sensitivity 
of the previous observation. 

4.2. The Molecular Outflow is still Moving Outward? 

In order to figure out the evolution of the molecular 
outflow OF, we compare the pressure between the dif- 
fuse molecular gas OF and the plasma gas in the CORE 
region. 

In the molecular outflow OF, the thermal pressure of 
constituent molecular cloud can be derived by P f,thm = 
n {kT, where n Q f is the number density of OF. Since we 
do not have density and temperature information of OF, 
we simply adopt typical values (100 - 1000 cm -3 and 10 

- 100 K) for CO(1-0) emitting gas. Thus the thermal 
pressure of OF is 1.4 x I0~( n ~ 13 ' dyne cm' 2 . 

In CORE region, the thermal pressures of the plasma 
outflow in the lower and the higher temperatures 
are lAf~ 1/2 x 10~ n dyne cm" 2 and 4.8/~ 1/2 x 
10' 10 dyne cm -2 , respectively (see Sec. 13.2. 3l and Tab. [2]). 
Adopting a volume filling factor of 0.1 - 0.01, their ther- 
mal pressures are (0.5 - 1.4) x 10~ 10 dyne cm -2 and (1.5 

- 4.8) x 10~ 9 dyne cm' 2 , respectively. Besides, the ram 
pressure that OF suffered from the plasma outflow can be 
derived by P ram = p p \ m x u 2 cl , where v rc i = |w oxp - u p i m 
is the relative velocity between the expansion velocity of 
OF, w cxp , and of the plasma outflow, v p i m . We assume 
that the expansion velocity of the plasma outflow is simi- 
lar to its thermal velocity, i.e., v p i m ~ y/ikT j /im p . Thus 
v p \ m estimated from the lower and the higher tempera- 
ture components in the CORE are 520 and 2,050 km s~ x , 
respectively. Since u oxp ~ 90 km s , v re \ for the lower 



13 



E 



o 
o 

_J. 

LU 

> 



200 



100 



-100 



-200 




-60 

POSITION OFFSET [arcsec] 



20 -40 -60 



Figure 12. Our NMA CO(l-O) p — v diagrams marked with the position of superbubbles claimed by Irwin & Sofuc ( 1996). The red curve, 
solid lines, and ellipse ar e the same with those in Figure[5] The big blue dot in each p — v diagram indicates the central velocity (v c ) of the 
superbubble claimed by Irwin & Sofuc (1996). The blue bar indicates the expansion velocity (v exp ) of superbubbles. The contour levels 
are 2, 3, 5, 10, 15, 20, and 50<r, where la = 9.09 mjy beam -1 , (a) The major-axis p — v diagram. Two superbubble A and C have the 
same central velocity of 55 km s — 1 . Their expansion velocities are 20 km s _1 and 150 km s _1 , respectively, (b) The p — v diagram parallel 
to the major axis with offset of —1.0". The superbubble B has a central velocity of —62 km s _1 and an expansion velocity of 33 km s — 1 . 
(c) The p — v diagram parallel to the major axis with offset of —5.0". The superbubble D has a central velocity of —126 km s -1 and an 
expansion velocity of 85 km s -1 . 
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Figure 13. The integrated intensity maps. The crosses mark the positions of superbubbles A, B, C, and D. (a) Integrated velocity range 
between 35 and 75 km s -1 . The contour levels are 2, 5, 10, 20, 40, 60, 80, 100, 150, and 200cr, where la = 219.45 mjy beam -1 km s -1 . 
(b) Integrated velocity range between 29 and 95 km s — . The contour levels are 2, 5, 10, 20, 40, 60, 80, 100, and 150(7, where la = 
194.75 mjy beam - 1 km s . (c) Integrated velocity range between —211 and —41 km s — . The contour levels are 2, 5, 10, 20, 40, 60, 80, 
and IOOct, where la = 229.70 mjy beam -1 km s — . 



and the higher temperature components in the CORE 
are 430 and 1,960 km s _1 , respectively, then the ram 
pressures of the lower and the higher temperature com- 
ponents are 2.9 J -1 / 2 x 10~ n dyne cm -2 and 1.3 / _1 / 2 
x 10~ 9 dyne cm -2 , respectively. After considering a fill- 
ing factor of 0.01 - 0.1, the ram pressures are (0.9 - 3.1) 
x 10~ 10 dyne cm -2 and (0.4 - 1.3) x 10~ 8 dyne cm -2 , 
respectively. We notice that in the lower temperature 
components, the ram pressure and the plasma thermal 
pressure have similar values, while in the higher temper- 
ature components, the ram pressure is higher than the 
plasma thermal pressure. Generally speaking, the ram 
pressures and the plasma thermal pressures in CORE 
are ~ 10 - ( 8 " 10 ) dyne cm -2 . Since the thermal pressure 
of the molecular outflow Of is ~ 10 - ( 11-13 ) dyne cm -2 
(Sec. 14. 2p . this indicates that the OF is pushed by the 
plasma gas. Thus we conclude that the OF is still mov- 
ing outward. 



4.3. The Evolutionary Stage of the Starburst Activity 

In order to figure out the evolutionary stage of the 
starburst activity, we defined an evolutionary parameter, 
r, which is the expansion timescale of molecular outflow, 



normalized to the total starburst timescale, fsB, 



to 



T = 



tsB ' 
t 



cxp 



(9) 
(10) 



where r = 1 means the starburst period is finished, and 
tsB includes two timescales, t exp , implying the timescale 
from the beginning of the starburst to now, and t cons , the 
molecular gas consumption timescale, i.e., the timescale 
that the current molecular gas in the starburst region 
to be used up, implying the timescale from now to the 
end of starburst period. t cons is highly depending on two 
processes, mass loss through molecular outflow, and star 
formation, 



SBR 



SFR 



(11) 



where Msbr is the molecular gas mass in the starburst 
region, and the SFR is the star formation rate. 

From Sec. 13.31 we know the mass loss rate, M Q f = 4.1 - 
8.5 M yr -1 . We define the size of the starburst region 
as the same with the base size of the X-ray emission, 
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Figure 14. The minor-axis p — v diagram of the superbubble D. 
This is averaged between —18" and —22" along m inor axis after 
rotati ng 17°, the same rotating degree with that of llrwin fe Sofuel 
(1996). The contour levels are 3, 5, 7, 10, 15, and 20cr, where la 
= 5.78 mjy beam -1 . The thin cross and the thin dashed circle 
indicate the central position/ velocity and t he expansion velocity , 
respectively, of superbubble D identified by Irwin & Sofuc (1996). 
From our NMA data, we re-identify the superbubble as D'. The 
thick dashed ellipse indicates the shell structure of superbubble 
D'. The thick cross indicates the central position of superbubble 
D'. 

which is ~ 0.5 kpc. Since the CO flux within starburst 
region is ~ 450 Jy km s _1 , this leads the molecular gas 
mass to be ~ 2.0 x 10 8 M Q . Since our data do not provide 
any information of the SFR, we use IR luminosity, Ltr, 
to estimate the SFR as the upper limit. We assume that 
the IR luminosity is all contributed from the central star- 
burst region, thus the SFR in the starb urst region can b e 
estimated as, SFR = 4.5 x 10~ 44 L m ([Kennicutd[l998l ). 
Taking L m = 10 1025 L (jSanders et all 120031 ) . SFR is 

3.2 M© yr" 1 . Thus t cons of NGC 3628 is calculated as 
17 - 27 Myr. 

Adopting a molecular outflow expansion timescale of 

3.3 - 6.8 Myr (from Sec. EES]), we have t SB = 20-34 Myr, 
and r = 0.11-0.25. This suggests that the evolutionary 
stage of the starburst activity is relatively young. 

5. SUMMARY 

Our NMA CO(1-0) observations detected a sub-kpc ~ 
370 - 450 pc scale molecular outflow for the first time 
with an expansion speed of ^ 90 ± 10 km s -1 in the 
starburst galaxy NGC 3628. The molecular outflow mass 
is 2.8 x 10 7 M Q , the expansion timescale is 3.3 - 6.8 Myr, 
and the mechanical energy is (1.8 - 2.8) x 10 54 erg. 

In order to understand the evolution of molecular 
outflow and the relation with the plasma outflow, we 
compared our NMA CO (1-0) data with the CXO X- 
ray archival data. We estimated the total energy both 
from mechanical energy of molecular outflow and ther- 
mal energy of plasma outflow. The total energy gen- 
erated by the starburst activity is estimated as (2.3 - 
2.8) x 10 55 erg. The thermal pressure of the molecular 
outflow is ~ io~( n ~ 13 ) dyne cm -2 . The thermal pres- 



sure and the ram pressure of the plasma outflow are ~ 
10~( 8 ~ 10 ) dyne cm~ 2 . This indicates that the molec- 
ular outflow is still pushed by the plasma gas and is 
moving outward. Besides, we confirmed one superbub- 
ble as a real s t ructu re from the four superbubbles that 
llrwin fe Sofuel (|1996l ) claimed. Two of them should be 
just noise, and the other one is difficult to be concluded. 

We also estimated the molecular gas consumption 
timescale, 17 - 27 Myr. This yields a total starburst 
timescale to be 20 - 34 Myr with an evolutionary param- 
eter of 0.11 - 0.25, suggesting that the starburst activity 
in NGC 3628 is still in a young stage. 
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Table 1 

Spectral Models for the CXO Spectra 



Region 
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13.5t^ a 


U - S1 -0.49 
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Note. — (1) Foreground absorbing column density of the VMEKAL in units of 10 22 cm 2 . (2) Temperature of component 1 in units of keV, 
where 1 represents the lower temperature component. (3) Normalization of component 1 in units of 10" 5 J 10~ 14 n c nudV f 4ir D 2 , where n c is 
the electron number density, rin is the hydrogen number density, V is the volume of plasma, and D is the distance between the sources and the 
observer. (4) Temperature of component 2 in units of keV, where 2 represent the higher temperature component. (5) Normalization of component 
2, the units is the same with that of component 1. (6) O abundance in units of Zq. (7) Mg abundance in units of Zq. (8) Si abundance in units 
of Zq. (9) Fe abundance in units of Zq. (10) Reduced x 2 of fit. (11) Degrees of freedom in the fit. 
a Fitted spectrum in the energy range 0.4 - 3.0 keV. 
b Fitted spectrum in the energy range 0.3 - 7.0 keV. 

Denotes parameters held fixed. 



Table 2 

The properties of plasma outflow in NGC 3628 



Region R Component kT EI n c M p i m Ppim.thm ^ P im,thm 
(1) (2) (3) (4) (5) (6) (7) (8) 

HALO 6.72 1 0.21 6.07 0.13 J" 1 / 2 40.0 J 1 / 2 0.09 J" 1 / 2 4.80 J 1 / 2 

2 0.60 3.30 0.09 f~ 1/2 29.4 J 1 / 2 0.18 f~ 1/2 10-1 f 1/2 

CORE 067 1 047 034 0.95 f~ 1/2 0.30 f i/2 1.43 f~ 1/2 0.08 f 1/2 

2 7.29 1.55 2.04 f- 1 / 2 0.64 f 1 / 2 47.7 f- 1 / 2 2.67 f 1 / 2 

Note. — (1) The radius of outflow in units of kpc. (2) The temperature component in absorbed 
two-temperature model, WABS* (VMEKAL + VMEKAL), where 1 represents the lower temperature 
component, and 2 represent the higher temperature component. (3) Temperature in units of keV. 
(4) Emission Integral (EI — f n e riHdV) in units of 10 61 cm" 3 . (5) Number density in units of 
10~ 2 cm -3 , where / is the filling factor. (6) Mass in units of 10 6 Mq. (7) Thermal pressure in units 
of 10" 11 dyne cm . (8) Thermal energy in units of 10 erg. 
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